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We present a numerical study of the density of states in a ferromagnet/superconductor junction and the Joseph- 
son current in a superconductor/ferromagnet/superconductor junction in the diffusive Umit by solving the Usadel 
equation with Nazarov's boundary conditions. Our calculations are valid for an arbitrary interface transparency 
and arbitrary spin-dependent scattering rate, which allows us to explore the entire proximity-effect regime. We 
first investigate how the proximity-induced anomalous Green's function affects the density of states in the fer- 
romagnet for three magnitudes of the exchange field h compared to the superconducting gap A: i) h^A, ii) 
/i > A, Hi) h ^ A. In each case, we consider the effect of the barrier transparency and allow for various con- 
centrations of magnetic impurities. We clarify features that may be expected in the various parameter regimes 
accessible for the ferromagnetic film, with regard to thickness and exchange field. In particular, we address how 
the zero-energy peak and minigap observed in experiments may be understood in terms of the interplay between 
the singlet and triplet anomalous Green's function and their dependence on the concentration of magnetic im- 
purities. Our results should serve as a useful tool for quantitative analysis of experimental data. We also investi- 
gate the role of the barrier transparency and spin-flip scattering in a superconductor/ferromagnet/superconductor 
junction. We suggest that such diffusive Josephson junctions with large residual values of the supercurrent at 
the O-TT transition, where the first harmonic term in the current vanishes, may be used as efficient supercurrent- 
switching devices. We numerically solve for the Josephson current in such a junction to clarify to what ex- 
tent this idea may be realized in an experimental setup. It is also found that uniaxial spin-flip scattering has 
very different effect on the O-tt transition points depending on whether one regards the width- or temperature- 
dependence of the current. Our theory takes into account vital elements that are necessary to obtain quantitative 
predictions of the supercurrent in such junctions. 

PACS numbers: 74.20.Rp, 74.50.+r,74.70.Kn 



I. INTRODUCTION 



Proximity structures consisting of ferromagnetic and su- 
perconducting materials are nowadays a very active research 
field. The interest in this type of systems has grown consider- 
ably over the last decade, since they offer novel and interesting 
phenomena to explore from a fundamental physics point of 
view. In addition, it is hoped that future applications in low- 
temperature nanotechnology may emerge from this research 
field. Ferromagnetism is usually considered to be antagonistic 
to conventional superconductors, since the exchange field acts 
as a depairing agent for spin- singlet Cooper pairs. However, 
closer examination reveals that the physical situation is more 
subtle than that. The proximity effect on a superconductor 
from a ferromagnet does not merely suppress the spin-singlet 
superconducting order parameter, but may also under specific 
circumstances induce exotic features such as odd-frequency 
pairing and long-ranged spin-triplet correlations^'^. 

Various theoretical idealizations allow for a relatively 
simple approach to ferromagnet/superconductor (F/S) het- 
erostructures in the quasiclassical framework. One of the most 
popular approaches in the literature employs the linearized 
Usadel^ equations with the Kupriyanov-Lukichev boundary 
conditions^, which is a viable method in the case of a weak 
proximity effect. This is for instance permissable when the 
barrier transparency of the F/S interface is low. Although the 
linearized treatment clearly represents a special limit, much 
useful information has been obtained through this approach. 



Among these are the oscillations of the anomalous Green's 
function^'^ in the ferromagnet due to the fact that the Cooper 
pair in the ferromagnet acquires a finite center-of-mass mo- 
mentum q = where h is the magnetic exchange en- 
ergy and is the Fermi velocity.^ Although several other 
works have also considered various aspects of the density of 
states in both diffusiveMdMLi^^min ^nd clean^^'^^'^^ ferro- 
magnet/superconductor junctions, most of these rely on sim- 
plifying assumptions concerning the interface, that it is either 
perfectly transparent or strongly insulating. 

From an experimental point of view, there has been sev- 
eral investigations of how the spin- splitting in energy level 
in a ferromagnet affects the proximity effect when placed in 
contact with superconductor^^'^^'^^'^^. Very recently, SanGior- 
gio et alP reported an anomalous double peak structure in 
the density of states of a Nb/Ni tunneling junction. They 
attempted to capture the qualitative features of the unusual 
subgap- structure of the density of states by using the Usadel 
equation, but were not able to do so. This clearly warrants fur- 
ther investigation, and serves as a motivation for employing 
more sophisticated models of a ferromagnet/superconductor 
interface, possibly including a domain structure in the ferro- 
magnet. 

Recently, a full numerical solution of the Usadel equation 
was employed by Gusakova et al including scattering on 
magnetic impurities, but under the simplifying assumption 
that the tunneling limit was reached at the F/S interface^"^. 
It would clearly also be of interest to consider higher bar- 
rier transparencies to elucidate how this may influence the 
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FIG. 1: (Color online). The ferromagnet/superconductor heterostruc- 
tures investigated in this paper. In a), a ferromagnet/superconductor 
bilayer setup is shown. We will model the barrier region at the 
normal/ferromagnet interface to be very strong, mimicking a free 
edge boundary condition or insulating region. Our proposed exper- 
imental setup is thus equivalent to the ones employed in the exper- 
imental works of Refs. [2Q||23| . In b), we consider a superconduc- 
tor/ferromagnet/superconductor junction. The superconductors are 
treated as reservoirs, thus unaffected by the proximity effect. We will 
consider a ferromagnetic layer thickness ranging from d/^ = 0.1 to 
d/^ — 1.0, where ^ is the superconducting coherence length. 



proximity effect. In this case, one should employ Nazarov's 
boundary conditions'^ instead of the Kupriyanov-Lukichev 
boundary conditions'^. Yokoyama et ah recently adopted 
this approach for the case without any magnetic impurities 
present'^'^^. When comparing the theoretical models against 
quantitative aspects of the experimental data, the influence of 
magnetic impurities and barrier transparency clearly play a 
pivotal role. For instance, as argued in Ref. 28, the predicted 
amplitude of the critical current in diffusive S/F/S junction 
is ~ 10^ times larger than the actual measured curves. This 
may be attributed to the strong suppression of Andreev bound 
states due to spin-flip scattering processes, which are usually 
not taken into account in the theoretical treatment (see, how- 
ever, Ref. Ilsl). To date, there exists no study of the density 
of states in F/S structures that allows access to both the full 
range of barrier transparencies and concentration of magnetic 
impurities. 

In the present paper, we remedy this by employing a nu- 
merical solution of the Usadel equation in the ferromagnetic 
region, taking into account an arbitrary concentration of mag- 
netic impurities as well as arbitrary barrier transparency. This 
permits us to comprehensively study the local density of states 
(DOS) in a dirty ferromagnet/superconductor junction for ex- 
perimentally realistic parameters. We study three cases for the 
size of the exchange field h compared to the superconducting 
gap A: i) h^A, ii) h > A, Hi) h:^A.ln each case, we com- 
pare the short-junction regime to the wide-junction regime 
to obtain the energy-resolved DOS. Among other things, we 



study the emergence of a zero-energy peak in the spectrum 
in a certain parameter range in addition to the transition from 
minigaps to peaks in the spectrum for increasing exchange 
field in the wide-junction regime. In particular, we investi- 
gate the interplay between the proximity-induced singlet and 
triplet anomalous Green's function in the ferromagnet and 
how these are affected by spin-dependent scattering. Since 
our results take into account an arbitrary proximity effect and 
magnetic impurity concentration, they should serve as a use- 
ful tool for performing quantitative analysis of experimental 
data. We envision a ferromagnet/superconductor bilayer as 
shown in Fig. [H which is virtually identical to the experimen- 
tal setup used in Refs. [2Qll23[ 

Another interesting issue in the context of ferromagnet- 
superconductor is that there have been put forth suggestions 
of exploiting superconductor/ferromagnet/superconductor 
(S/F/S) Josephson junctions as qubits in quantum computers. 
By now, the O-tt oscillations that occur in S/F/S junctions are 
also wefl-estabfished both theoretically^^'^^'^''^^'^"^'^^'^^ and 
experimentallyi^^'Q^^^^^^^^^^^Q^^i. 

As the fundamental understanding of the physics in an 
S/F/S junction begins to take shape, the next step to take 
should be a more sophisticated modelling of such systems in 
order to achieve better quantitative agreement between theory 
and experiment. Parameters such as barrier transparency and 
magnetic impurities become important in this respect, as they 
can have crucial impact on the behaviour of the supercurrent 
in an S/F/S system. To this end, several recent works inves- 
tigated some aspects of the supercurrent by including mag- 
netic impurities'^'^^'^^'^^'^^'^^'^^, but restricting themselves to 
the tunneling limit. Under the opposite simplifying assump- 
tion of transparent interfaces, the influence of spin-flip scat- 
tering has also been investigated"^^. The Josephson current in 
an S/F/S junction may under these circumstances be written 
to good accuracy as / = /q sin (j), where (j) is the macroscopic 
phase difference. As a direct result, the usual O-tt oscillations 
observed in S/F/S junctions do not exhibit any residual value 
of the supercurrent at the transition points where the current 
switches sign. 

However, an interesting opportunity presents itself in the 
case when the current-phase relationship deviates strongly 
from its usual sinusoidal behavior. Higher harmonics in the 
Josephson current will lead to a finite residual value right at 
the cusps of the critical current oscillations with exchange 
field or temperature, and these cusps are indicative of a sign 
reversal of the supercurrent. This simple observation gives 
rise to a very interesting prospect. If the residual value of 
the current at the cusps were to be of considerable magnitude, 
one could exploit this effect to create a dissipationless current 
switching device. In this case, the direction of the Joseph- 
son current could be instantaneously flipped by some external 
control parameter when the system parameters are such that 
the junction is very close to the O-tt transition point. 

Such an idea depends, however, on the possibility to actu- 
ally obtain a large residual value of the current near the cusps. 
The early experiments ^^'^^'^^'^^^^^ measuring the Josephson 
current in diffusive S/F/S systems reported that the critical 
current vanished at these cusps. But recently, Sellier et al^^- 
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observed a finite albeit small supercurrent at the O-tt transition 
point where the first harmonical term vanishes. 

In Ref. 28, the magnitude of the residual value of the cur- 
rent at the cusps was investigated by solving the linearized Us- 
adel equation^ using Nazarov's boundary conditions'^. How- 
ever, this calculation was performed assuming two limiting 
circumstances: i) no magnetic impurities and ii) a weak prox- 
imity effect. The residual value was also analyzed in the clean 
limit in Ref. It is definitely of much interest to go beyond 
these approximations and see how the physical picture is al- 
tered. As mentioned previously, a consequence of i) is that 
the theoretically predicted magnitude of the critical current is 
a factor 10^ larger than the experimentally measured value. 
Even though highly valuable qualitative information may be 
obtained in approximations such as i) and ii), the progress on 
both the theoretical and experimental side of S/F/S Josephson 
junctions calls for a higher accuracy of the quantitative predic- 
tions. Obviously, this is also of paramount importance in the 
context of discussing a practical supercurrent- switch device. 

We here develop a theory for the supercurrent in an S/F/S 
junction which takes into account an arbitrary concentration 
of magnetic impurities and arbitrary transparency of the in- 
terfaces. Both of these are of vital importance in obtaining a 
quantitative agreement with experimental findings. We nu- 
merically solve the Usadel equation and employ Nazarov's 
boundary condition. In the intermediate barrier transparency 
regime, we find that a finite but small residual value of the 
current is permitted at the O-tt transition points. However, the 
effect of spin-flip scattering is very different when comparing 
the width-dependence against the temperature-dependence of 
the current. The transition point is much more robust in the 
width-dependence, although increasing the concentration of 
magnetic impurities reduces the residual value. On the other 
hand, increasing spin-flip scattering completely removes the 
O-TT transition point for the temperature-dependence. 

This paper is organized as follows. In Sec. [Ill we es- 
tablish the theoretical framework which will serve as our 
tool to obtain both the density of states in the diffusive fer- 
romagnet/superconductor bilayer and the Josephson current 
in the diffusive superconductor/ferromagnet/superconductor 
junction. In Sec. Hill we present our results with a discus- 
sion of these in Sees. IIII ASI and llll A6[ as well as in Sec. |lVl 
Finally, we summarize the results in Sec. |Vl 



II. THEORETICAL FORMULATION 

The central quantity in the quasiclassical theory of 
superconductivity is the quasiclassical Green's functions 
^(pF,R;£:,t), which depends on the momentum at Fermi 
level pf, the spatial coordinate R, energy measured from 
the chemical potential e, and time t. A considerable lit- 
erature covers the Keldysh formalism and non-equilibrium 
Green's functions^^*^^-^^*^^*^^. Here we only briefly sketch 
the theoretical structure, for the sake of readability and for 
establishing notation. The quasiclassical Green's functions 
^(Pf, R; t) is obtained from the Gor'kov Green's functions 
G(p, R; £, by integrating out the dependence on kinetic en- 



ergy, assuming that G is strongly peaked at Fermi level, 

^(PF, R; t) = ^ J depG(p, R; e, t). (1) 

The above assumption is typically applicable to superconduct- 
ing systems where the characteristic length scale of the per- 
turbations present, namely superconducting coherence length, 
is much larger than the Fermi wavelength. The correspond- 
ing characteristic energies of such phenomena must be much 
smaller than the Fermi energy £f. The quasiclassical Green's 
functions may be divided into an advanced (A), retarded (R), 
and Keldysh (K) component, each of which has a 4 x 4 matrix 
structure in the combined particle-hole and spin space. One 
has that 

where the elements of ^(pf, R; t) read 

The quantities g and / are 2 x 2 spin matrices, with the struc- 
ture 

- \9i] 9ii) 

Due to internal symmetry relations between these Green's 
functions, all of these quantities are not independent. In par- 
ticular, the tilde-operation is defined as 

/(pF,R;^,t) = /(-pF,R;-e,t)*. (5) 

The quasiclassical Green's functions ^(pf, R; ^, t) may be de- 
termined by solving the Eilenberger^^ equation 

[eps -t,g]^^ ivfV^ = 0, (6) 

where S contains the self-energies in the system such as im- 
purity scattering, superconducting order parameter, and ex- 
change fields. The star-product is noncommutative and 
is defined in Appendix. When there is no explicit time- 
dependence in the problem, the star-product reduces to normal 
multiplication. This is the case we will consider throughout 
the paper. The operation p^^g inside the commutator should 
be understood as = diagjpa, pa}^. Pauli-matrices in 
particle-hole X spin (Nambu) space are denoted as pi, while 
Pauli-matrices in spin-space are written as r^, all of which are 
defined in the Appendix. The Green's functions also satisfy 
the normalization condition 

g®g = l. (7) 

In the special case of an equilibrium situation, one may ex- 
press the Keldysh component in terms of the retarded and ad- 
vanced Green's function by means of the relation 

^^ = (^^-^^)tanh(/3£/2), (8) 
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where f3 = is inverse temperature. In nonequilibrium sit- 
uations, one must derive kinetic equations for nonequilbrium 
distribution functions in order to specify the Keldysh part^^. 

The above equations suffice to completely describe for in- 
stance a single superconducting structure, but must be sup- 
plemented by boundary conditions when treating heterostruc- 
tures such as F/S junctions. These boundary conditions take 
different forms depending on the physical properties of the in- 
terface. The Kupriyanov-Lukichev^ boundary conditions may 
be applied for a dirty junction in the tunneling limit when the 
transparency of the interface is low, while an arbitrary inter- 
face transparency requires usage of the boundary conditions 
developed by Nazarov^^. Boundary conditions for a spin- 
active interface have also been derived^^-^^. 

We will consider the dirty limit of the Eilenberger equation 
Eq. S, which leads to the Usadel equation^. Our motiva- 
tion is that this is the experimentally most relevant situation. 
This will be an appropriate starting point for diffusive systems 
where the scattering time due to impurities satisfies Xr <C 1, 
where X is the energy scale of any other self-energy in the 
problem. For strong ferromagnets where h becomes compa- 
rable to ^F, the stated inequality may strictly speaking not 
be valid for X = h. Hence, we will restrict ourselves to 
the regime h/e-p <C 1 and assume that^^ /ir <C 1. Below, 
we will mostly concern ourselves with the retarded part of 
^(Pf, R; £, t), since the advanced component may be found 
via the relation 



-(P3/P3)^ 



(9) 



By isotropizing the Green's function due to the assumed fre- 
quent impurity scattering, it is rendered independent of pf. 
This isotropic (in momentum space) Green's function satis- 
fies the Usadel equation in the ferromagnet: 

DV{gVg) + i[ep^ M - a,f - ^so, g] = 0. (10) 

Above, the exchange energy h is accounted for by the ma- 
trix M = dmg{hTs, hrs), assuming a magnetization in the 
z-direction, while the spin-flip self-energy reads 

oTsf ^ 

(jso(R;e) = -^^aiP3g{'R]e)p3ai, (11) 

where r^f is the spin-flip scattering time and Si is a spin ex- 
pectation value, while r^o is the spin-orbit scattering time. We 
have defined the matrices ai = diag(Tj_, t^^). The diffusion 
constant is given hy D = Vp'^imp/S, where rimp is the impurity 
scattering relaxation time. We will here consider either uniax- 
ial spin-flip scattering, such that ^'3 = 1 (z-direction) and zero 
otherwise, and also isotropic spin-flip scattering where Si = 1 
fori e {1,2,3}. For later use, we denote Sxy = Si = S2 and 
S, = Ss. 

Let us now consider the retarded part of Eq. ([TOb which has 
the same form, namely 

DVig'^Wg^) + i[ef3 + M - a^f - g""] = 0. (12) 



From now on, we will omit the superscript 'R' on the Green's 
function. Moreover, we will find it convenient to parametrize 
the Green's function by exploiting the normalization condi- 
tion. In the 5-wave superconductor {x > d) and normal metal 
(x < 0), we use the bulk solutions 



1 Q 

-1 



(13) 



where g = cosh(6>s), / = sinh(6>s), Os = atanh(A/£:). 
The Green's function in the ferromagnet may conveniently be 
parametrized as^ 



9f 



I cosh6>|(£:) 



V— sinh^|(£) 











cosh^l(6:) 
-sinh l9|(£) 




sinh(9|(£:) 
- cosh 6^1 (e) 




sinh^l(£:) 



- COSh^|(£)y 

(14) 

which satisfies g\ — \. We have made use of the symme- 
try d^{e) = 6>*(-£) in obtaining Eq. ([T4|). Note that for 
= 0, d^{e) = -Oi{s) is satisfied. Note that in Eq. ([T4l) . 
there is both a singlet and triplet component of the anoma- 
lous Green's function. The triplet component is opposite- 
spin paired {Sz = 0), and there are no equal-spin pairing 
{Sz = ±1) components in our system since we consider 
homogeneous magnetization and a spin-inactive barrier. The 
Sz = triplet component nevertheless plays a pivotal role in 
interpreting the behaviour of the density of states, as we shall 
see later. This is because it has a special symmetry property 
refered to as odd-in-frequency, which will be elaborated upon 
in Sec, nil A 41 

The Usadel equation Eq. ([T2l) then yields for majority and 
minority spin (a =t, i= ±1) in the ferromagnet: 



DdlO^ ^ 2i{e ^ ah) sinhOa 



crSx 



2rsi 



■sinh(^^-^;) 



S 1 

- sinh2i9^ sinh((9t + 6>|) = 0. (15) 

4rsf 2rso 

We now use the generalized Nazarov boundary condition valid 
for the diffusive regime^^. At x = {0, (i}, it reads 



'^iL.RLgpdxgF 



= T 



4 + r({^F,^L,i?}-2) 



= {0,d} ^ ' ■ '^^'^-^J / x={0,d} 

(16) 

where [. . .] and {. . .} denote the commutator and anticom- 
mutator, respectively, and gL,R{= gN,s) denotes the Green's 
function on the left and right side of the ferromagnet. The 
parameter "Jl^r = Rb'^/Rf denotes the ratio between the 

L R 

resistance in the left/right barrier region R^' and the resis- 
tance in the F region Rp. We have conventionally introduced 
the parameter r, which is the transmissivity of the interface^^. 
Giving an expression for r in terms of microscopic parameters 
of the interface is not very practical, and we will therefore use 
r as a phenomenological parameter to characterize the trans- 
parency of the interface, in complete analogy with Eq. 36 of 
Ref. I25i. Here, r = corresponds to zero transmission of 
quasiparticles incident on the superconducting interface, and 
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r = 1 corresponds to perfect transmission. Below, we will 
consider two values of r, namely 0.1 corresponding to low 
transmissivity, and r = 0.5 corresponding to intermediate 
transmissivity. Note that the parameters r and 7 may be varied 
independently^^, and one does not in general have 7 ~ r~^. 
The reason for this is that 7 is related to the constriction area 
of the junction, which may be altered independently of the 
scattering strength proportional to r~^. Inserting Eqs. ([T3l) 
and ([T4b into Eq. ([T6b yields the boundary conditions 



2r sinh6> 



2 — r + r cosh 0^ ' 

2r sinh^i 
2 — r + r cosh 6^ ' 



(17) 



at the normal/ferromagnet (N/F) interface {x = 0), while at 
the F/S interface {x = d) we have 



InddxO^ - 

Z — T 



2r(cosh O^f — sinh O^g) 
- r + r (cosh 6^g — sinh /) ' 
— 2r(cosh^l/ + sinh^i^) 
2 - r + r(cosh Oig + sinh OJ) ' 



(18) 



For later use, we define the Thouless energy et = D/d^. As 
a measure of the strength of the spin-flip and spin-orbit scat- 
tering, which increases with decreasing spin relaxation time, 
we introduce ^sf = 'tT^ and ^§0 = • Also note that the 



1 for isotropic spin-flip scattering, while 5'^ = 3 



and Sxy = for uniaxial spin-flip scattering along the z- 
direction. The spin-resolved and normalized DOS is obtained 
as 

= Re{cosh i9|}, = Re{cosh(9|}. (19) 

Furthermore, we define the total DOS sls N = ^^N(j/2. 
Eqs. ([T5]) , (fTTl) , and (fTSl) now constitute two coupled non- 
linear second order differential equations supplemented with 
boundary conditions which may be solved numerically. 

Before presenting our results for the DOS, we establish 
the theoretical framework for our treatment of the Joseph- 
son current in an S/F/S junction with spin-dependent scatter- 
ing, following the notation of Ref. 61. The physical system 
studied here consists of a junction with two 5-wave super- 
conductors separated by a diffusive ferromagnet with a re- 
sistance Rf and length d much larger than the mean free 
path. We here only consider uniaxial spin-flip scattering. 
The interface regions are characterized by a resistance Rb- 
The transparencies of the junction interfaces are given by 
T = 4 cos^ V^/(4 cos^ + Z^) where Z is a measure of the 
barrier strength, and the barriers themselves are considered to 
be spin-inactive and modelled by infinitely narrow insulating 
barriers [/(x) = {Zvf /2)[S{x — d) -\- S{x)]. Above, 'i;^ is the 
Fermi velocity and (p is the injection angle measured from the 
interface normal to the junction. We employ the quasiclassi- 
cal theory of superconductivity^Si^iiS^iSiSiS^, and make use of 
the 6>-parametrization^^ of the Green's function. The retarded 
part g^ may, due to symmetry requirements, be written as^^ 



g = sin 6>(cos ipri + sin i/jt^) + cos 6^73^. 



From the Usadel quation, one obtains 

D[dlO - {^x^pf cos sin 0] + 2i(£: + a/i + 17) sin l9 = 0, 
dx[sm^O{dx^)]=0, (21) 

for a =T,i spins. Here, D is the diffusion constant, h is 
the exchange field, and 7 is the self-energy associated with 
uniaxial spin-flip scattering. We employ the bulk solutions 
of the Green's function in the superconducting regions, as- 
suming that these are much less disordered than the ferromag- 
netic layer. To gain access to the full regime of different bar- 
rier transparencies, we again make use of Nazarov's boundary 
conditions'^ which are valid for a non-magnetic but otherwise 
arbitrary contact. These boundary conditions at x = {0, 
may in the present case then be written as^^ 

L{RB/RF)dxO = ±/±[J^cosi9cos(V^ t ^/2) -GsinO], 
L{Rb/Rf) sinOdxtp = tI±^^H^ T ^/2), (22) 

where the upper (lower) sign is valid at x = d (x = 0). We 
have defined the following quantities: 

/± = (4T/[A±(1 + ^2 + a')]), A± = (2 - T) 

+ 2T [J^ sin l9 cos (^ T 0/2) + S cos6>]/(l + + S^). 

(23) 

and denoted the phase in the right and left superconductor as 
±0/2, giving rise to a total phase difference of 0. Above, 



(24) 



denote components of the bulk Green's function in the su- 
perconductor. Eqs. (I2TI) and (1221) constitute two coupled, 
second order, nonlinear differential equations with bound- 
ary conditions and may be used to numerically solve for 
= 0{x^e). Once 6{x^e) is obtained, the retarded Green's 
function is specified everywhere in the diffusive ferromagnet. 
The current-density may then be calculated by 



j = -{NF\e\D/4) / deTr{ps{gdxgr}, (25) 

^ — CO 

where Nf is the density of states at Fermi level in the nor- 
mal state, |e| is the electronic charge, ps is a Pauli ma- 
trix in particle-hole space, g is the full Green's function in 
Keldysh(g)particle-hole(8)spin space, and the superscript 'K' 
denotes the Keldysh component. The current is obtained by 
I = JSa, where is the surface area of the junction, and 
the critical current is defined as Ic = max{/(0)}. The IqRn- 
product can now be computed, where Rn = 2Rb -\- Rf- In 
the actual calculations and numerical implementation, we em- 
ploy the Matsubara representation e ^ ico, and parametrize 
the quasiclassical Green's functions with the quantity as 



fuj = guj^uj/^ sinl^ 



-i'0 



r-^ = ^l^gju; = sin Oe'^, 



cos^. 



(26) 



(20) 



We may now express Eqs. (I2TI) and (1221) in terms of the 
fermionic frequency co = (2n + l)7r//3, n = 0, 1, 2 . . ., and 
the Green's functions in Eq. 
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For later purposes, we introduce dimensionless measures of 
the Thouless energy and inverse spin-flip scattering lifetime: 
£ = StI ^ = D/{(PA) and g = 7/ A. For simplicitly, we 
will also neglect the spatial depletion of the superconducting 
order parameter near the interfaces. 

III. RESULTS 

We now proceed to present our results for the DOS and 
Josephson current in detail. In the first part of this section, 
we consider a diffusive ferromagnet/superconductor bilayer, 
while in the second part we investigate a diffusive supercon- 
ductor/f erromagnet/superconductor j unction . 

A. Density of states in ferromagnet/superconductor junction 

We will divide our results for the DOS into subsections to 
clarify the role of the spin-dependent scattering for different 
exchange fields and film thicknesses. In the spirit of Ref. llTL 
we will consider the cases i) h^A, ii) h > A, in) A. 
For h A, the proximity effect in the ferromagnet is weak 
unless the Thouless energy is very high {st ^ A). Through- 
out the rest of this paper, we fix 71, = 100 and = 1. 
This corresponds to a scenario where the normal metal reser- 
voir effectively acts as a very strong insulating barrier, mim- 
icking a vacuum boundary. Since the resistance of the N/F 
interface will constitute the largest contribution to the total re- 
sistance, the conductance of the junction will be equivalent 
to the DOS at the N/F interface in the tunneling limit. We 



will study an intermediate value r = 0.5 of the barrier trans- 
parency, since this cannot be reached within the usual approx- 
imations of either a fully transparent interface (r = 1) or a 
tunneling barrier (r <C 1), and contrast this with a strong 
barrier r = 0.1. These choices for the barrier strength are 
experimentally the most relevant ones. For each of the cases 
i)-in), we will investigate two different thicknesses of the fer- 
romagnetic layer, namely d/(_ = {0.1, 1.0}, corresponding to 
Et/A = {100, 1}. Here, ^ = y^D/A is the superconduct- 
ing coherence length. We comment further on our choice of 
parameters in Sec. [iVl For each type of spin-dependent scat- 
tering, we define the dimensionless parameter as a measure 
of the inverse scattering time. For spin-flip scattering we have 
^so = and g = g^i/sT with Sz = 3, Sxy = in Eq. ([T5l) for 
the uniaxial case while Sz = Sxy = 1 in the isotropic case. 
For spin-orbit scattering, we have ^sf = and g = gso/^T- 
Unless otherwise specified, the DOS is calculated right at the 
interface between the normal metal and the ferromagnet, i.e. 
X = (see Fig. [T]). This corresponds precisely to the experi- 
mental situation in Ref . [13. 

In what follows, we first present our numerical results for 
the DOS for the cases i)-m) described above. We then in- 
vestigate and explain the features seen in each of those cases 
in separate subsections. One of the main conclusions in this 
section is that the distinction between different types of spin- 
dependent scattering, e.g. spin-flip and spin-orbit scattering, 
may actually become very important in terms of interpreting 
the DOS in a ferromagnet/superconductor bilayer. We relate 
this to the behaviours of the proximity-induced anomalous 
singlet and triplet Green's functions in Sec. IIII A5[ 



1. Case i): h<A 

We first consider the case of a very weak exchange field, 
h/A = 0.5, which splits the DOS for spin-| and spin-| 
electrons. A peculiar feature that may arise in a ferromag- 
net/superconductor junction is an enhancement of the DOS at 
zero energy, manifested as a peak. This issue was addressed 
in Ref. .17.26^27,. and has also been experimentally observed 
in Ref. [2O. As seen in Fig. [21 the ferromagnet becomes 
fully proximized by the superconductor in the thin-layer case 
d/^ = 0.1, and the DOS is equivalent to the bulk of the su- 
perconductor (in the absence of spin-dependent scattering). In 
fact, the DOS is virtually unaltered compared to the paramag- 
netic case for the short-junction regime regardless of whether 
we consider the tunneling or intermediate transparency case. 
However, an important distinction between different types of 
spin-dependent scattering becomes evident in the thin-layer 
case d/^ = 0.1. The effect of spin-flip scattering, whether it 
is uniaxial or isotropic, is that the DOS-gap closes upon in- 
creasing g. However, increasing the spin-orbit scattering rate 
does not affect the DOS in any way. We investigate this issue 



in more detail in Sec. IIII A 61 

The situation changes dramatically when going to a wider 
junction regime, here modelled by d/^ — 1.0. In the tun- 
neling limit, a minigap-like structure opens up in the DOS 
at energies e = ±h. This can be understood as a result of 
the exchange- splitting in the ferromagnet which shifts the en- 
ergies, and hence the density of states, for majority and mi- 
nority spins hy ±h. Increasing the barrier transparency to 
r = 0.5, the minigap-like structure is retained and widened, 
which is reasonable since the proximity effect becomes larger 
and the minigap scales with r. Another feature that is seen in 
the d/^ = 1.0 case is the appearance of a zero-energy peak 
(ZEP) in the DOS. The peak protrudes with increasing spin- 
flip scattering and is actually split into two in the absence of 
magnetic impurities. In Ref.^"^, where the tunneling limit was 
applied to the F/S interface, no zero-energy peak was found 
to appear in the spectra. In the clean limit, however, a zero- 
energy peak was found to appear in the DOS for a wide range 
of parameters'^. Although the appearance of a ZEP has been 
investigated in previous work, none of these considered the 
effect of spin-dependent scattering. In order to understand the 
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FIG. 2: (Color online). The local density of states at x = (see Fig. [T]for /i/A = 0.5. For each thickness of the ferromagnetic layer and 
barrier transparency, we study the role of different types of spin-dependent scattering. 



ZEP feature seen in Fig. [21 we study the dependence of the 
anomalous Green's functions on spin- selective scattering in 
Sec, nil A 51 



DOS in the ferromagnetic layer to deviate less from unity. 
3. Case in): h ^ A 



2. Case ii): h> A 

Next, we consider the case h/ A = 1.5, shown in Fig. [S] 
As in the previous case, the DOS is almost unaffected by 
the presence of an exchange field in the short-junction regime 
d/^, = 0.1. Also, the gap in the DOS shows a remarkable re- 
siliance towards increasing the spin-orbit scattering. However, 
for d/^ = 1.0 the DOS is strongly modified and displays two 
peaks located ate = ±A and e = ±h, respectively. Increas- 
ing the concentration of magnetic impurities suppresses these 
peaks. In the clean limit, Zareyan et al found a similar de- 
velopment of the DOS with increasing Thouless energies (see 
Fig. 4 in Ref.[l7|). Note that the suppression of the proximity- 
induced features in the DOS due to the superconductor is now 
stronger for a given (i/^ as compared to the case h < A (Fig 
O. In general, increasing the value of h leads to a smaller 
magnitude of the proximity-induced anomalous Green's func- 
tions (see also Fig. [7]), which in turn causes the normalized 



Finally, we investigate the case h A. If the ferro- 
magnetic layer is an alloy of the type Cui-a^Ni^^, a reason- 
able value of the exchange field may be found in the range 
h = 10 — 50 meV. Here, we consider h/A = 15 and 
h/ A = 50. With increasing value of the exchange field, the 
proximity effect becomes weaker. Therefore, we restrict our 
attention to the short-junction regime (d/^ = 0.1) with an 
intermediate barrier- transparency (r = 0.5) since the DOS 
deviates little from unity in the tunneling regime and for wide 
ferromagnetic layers. 

For h/A = 15, the usual peak at e = A is present, and 
a minigap structure is seen in the absence of spin-flip scat- 
tering. Interestingly, increasing spin-flip scattering not only 
closes the minigap, but actually causes the DOS to develop a 
peak at zero energy. This may be understood by considering a 
subtle interplay between the singlet and triplet components of 
the proximity-induced anomalous Green's function in the fer- 
romagnetic layer, and we delay a detailed explanation to Sec. 
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FIG. 3: (Color online). The local density of states at x = (see Fig. [T]for hj IS. — 1.5. For each thickness of the ferromagnetic layer and 
barrier transparency, we study the role of different types of spin-dependent scattering. 



IIIIA5I 

For /i/ A = 50, all peaks and minigap features are now ab- 
sent for I el < A, and the only feature remaining in the spec- 
trum is a dip at £ = A. In this case, the qualitative effect of 
the different types of spin-dependent scattering is the same. 
Upon increasing djS^ even further (d/^ ^ 1), correspond- 
ing to a weaker proximity effect, one finds that the correction 
to the DOS oscillates upon increasing the spin-flip scattering 
rate, in contrast to the monotonous decay that might have been 



expected^"^. These oscillations are a result of the modified os- 
cillation length of the proximity-induced anomalous Green's 
function in the F region 

Also note that although the corrections to the normal- 
state DOS diminish rapidly iox h ^ IS. upon increasing 
d/^, the combination of lock-in detection with an ultra-low 
noise DC/ AC mixer permits the resolution of structures in 
the DOS up to a factor 10~^ smaller than the background 
conductance^^. 



4. Zero-energy peak 



It is instructive to consider the development of the zero- 
energy peak in the DOS. In Fig. [5l we plot the DOS around 
zero energy for d/^ = 1.0 for increasing values of the ex- 
change field h. The gradual formation of a zero-energy peak 



is clearly observed, and the peak flattens out when h becomes 
sufficiently large. The physical reason for this is intriguing. 
Yokoyama et alP have related this phenomenon directly to 
the proximity-induced odd-frequency pairing component of 
the anomalous Green's function in the ferromagnet (see also 
Refs.^^'^^). In these works, it is shown that the DOS in the fer- 
romagnet is enhanced due to the emergence of the proximity- 
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FIG. 4: (Color online). Local density of states at x = (see Fig. [T]) for /i/A = 15 (upper row) and /i/A = 50 (lower row). Here, we have 
fixed d/^ = 0.1 and r = 0.5 and study the role of different types of spin-dependent scattering. 



induced odd-frequency pairing. Since we are considering 
the isotropic part (with respect to momentum) of the Green's 
function due to the angular averaging in the dirty limit, one 
would perhaps naively expect that only the singlet compo- 
nent should be present. This is because the singlet anomalous 
Green's function is usually taken to be even under inversion 
of momentum, while the triplet components are taken to be 
odd under inversion of momentum. However, another pos- 
sibility exists that allows for the presence of triplet correla- 
tions in the ferromagnet, involving a change of sign of the su- 
perconducting order parameter under inversion of frequency. 
This type of pairing has been dubbed odd-frequency pairing 
in the literature.^ Inversion of frequency is equivalent to an 
exchange of (relative) time coordinates for the field operators, 
since e is the Fourier transform of the relative time coordinate 
t = ti—t2- Note that although even-frequency triplet correla- 
tions are destroyed in the dirty limit due to the isotropization 
stemming from impurity scattering, odd-frequency triplet cor- 



relations may persist, since these do not vanish under angular 
averaging. 

It remains to be clarified how the ZEP in a ferromag- 
net/superconductor junction is affected by spin-dependent 
scattering. To investigate this, we plot the DOS in Fig. [6l 
for a fixed exchange field /i/A = 0.3 which gives a ZEP 
in the absence of spin-dependent scattering, and then succes- 
sively increase the scattering rate. It is seen that the effect 
of increasing the spin-flip scattering rate (both uniaxial and 
isotropic) is a suppression of the proximity-induced features 
in the DOS. Qualitatively, the same occurs upon increasing 
the spin-orbit scattering rate (right panel of Fig. [S]), but an in- 
teresting difference is that the peak is eventually transformed 
into a dip at e = 0. Increasing the spin-orbit scattering rate 
further {g ^ 1) leads to a fully developed minigap in the 
DOS. We propose an explanation of this peculiar phenomenon 
in the following section. 



5. Anomalous Green's functions 

In order to understand the interplay between the singlet and 
triplet components of the induced superconducting anomalous 
Green's function in the F region with regard to the zero-energy 
behaviour of the DOS, consider the anomalous Green's func- 
tions defined as: 

ft{e,x) = [sinh l9| (£, x) +sinhl9|(6:,x)]/2. (27) 

At zero energy e = 0, one finds that Re{/s}=lm{/t}=0. 
The reader is reminded of the relation between the anoma- 



lous Green's functions and the DOS, which is a physical ob- 
servable: the total DOS is given as = Re{cosh6>cr}/2. 
Also note that the singlet and triplet components differ not 
only in their spin- symmetry, but also with respect to their 
energy-dependence (even and odd, respectively), as noted in 
Sec. lIIIA4l 

In Fig. [71 we plot Imj/g} and Re{/t} as measures of the 
singlet and triplet induced proximity Green's function in the 
ferromagnet right at the N/F interface {x = 0). We consider 
the effect of each type of spin-dependent scattering separately. 
As seen, the peak in the singlet component occuring at a finite 
value of h vanishes upon increasing spin-flip scattering rate. 
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FIG. 5: (Color online). Plot of the DOS of a ferromag- 
net/superconductor junction with increasing exchange field h, using 
d/^ — 1.0 and r = 0.5. Here, we have set the spin-dependent 
scattering to zero. 

and the decay eventually becomes monotonuous for ^^0.3. 
However, this is not the case for the triplet component: the 
maximum value of |Re{/t}| occurs at a finite value of h even 
upon increasing the spin-dependent scattering rate. Therefore, 
the triplet component may become similar in magnitude to 
the singlet component even for weak exchange fields if the 
spin-flip scattering rate is sufficiently large. This explains the 
appearance of a zero-energy peak in the plot for d/^ = 1.0 in 
Fig. [3] upon increasing g. 

Another interesting feature can be observed in Fig. [71 From 
the plots, it is seen that the anomalous Green's functions (both 
singlet and triplet) for a given exchange field may actually 
become larger upon increasing the rate of spin-flip scatter- 
ing. This becomes evident around h/ /A = 0.1 when com- 
paring the lines for ^ = 0.0 with g = 0.1. Thus, we have at 
hand the opportunity to see an enhanced proximity-effect by 
spin-flip scattering due to the increased magnitude of \fs^t \ at 
low values of h. The enhanced proximity-effect actually be- 
comes very pronounced also in the case of bulk odd-frequency 
superconductors^"^ . 

From Fig. [71 we may actually also find an explanation 
for the remarkable behaviour of the DOS upon increasing the 
spin-orbit scattering rate in Fig. [6l For concreteness, let us 
first focus on the regime h/A 0.3 corresponding to Fig. 
[6l Let us compare the plots for spin-flip and spin-orbit scat- 
tering. For both uniaxial and isotropic spin-flip scattering, it 
is seen that the singlet component does not change much in 
magnitude, while the triplet component is suppressed. How- 
ever, the singlet and triplet components are still comparable 
in magnitude for ^ = 0.5, which accounts for the suppression 
of the minigap feature which is due to the singlet component. 
The situation is markedly different for the spin-orbit scatter- 
ing. Now, the singlet component is actually enhanced in mag- 
nitude (for /i/ A :^ 0.3) while the triplet component is reduced 
very strongly. In fact, the triplet component becomes essen- 
tially zero for very large spin-orbit scattering rates, while the 
singlet component still may have a considerable magnitude. 



This explains the appearance of a minigap in Fig. [6l upon 
increasing the spin-orbit scattering rate: the singlet compo- 
nent, which is responsible for the minigap, increases while 
the triplet component decreases. 

The above discussion was restricted to /i/ A :^ 0.3 for con- 
creteness, but the results nevertheless allude to a much more 
general principle: the spin-orbit scattering is much more detri- 
mental for the triplet component than the singlet component. 
To elucidate this feature, consider Fig. [8l where we com- 
pare the proximity-induced singlet and triplet component in 
the ferromagnetic layer. We set d/^ = 1.0, since for d/^ <C 1 
the singlet component completely dominates the triplet com- 
ponent (see Fig. [2land[3l). Also, we set r = 0.5, but we 
checked that the qualitative features are identical for lower 
barrier transparencies. As seen from the figure, the singlet 
component is much more robust towards increased spin-orbit 
scattering than the triplet component. The latter only becomes 
appreciable in magnitude when the exchange field becomes of 
the same order or larger than A when the spin-orbit scattering 
range is large, ^ 1. Note that in contrast, both the singlet 
and triplet components are strongly reduced with increasing 
spin-flip scattering. 

In the absence of an exchange field, it is clear from Fig. [3 
that the singlet component is completely independent of the 
rate of spin-orbit scattering. Although we have only shown 
this explicitly for £ = 0, we have numerically confirmed that 
the magnitude of the singlet component at = remains 
unchanged upon increasing g. On the other hand, the triplet 
component is always zero at h = 0, but remains very close 
to zero upon increasing the exchange field h for high values 
of ^. A natural question arises: why is the singlet component 
insensitive to spin-orbit scattering while the triplet component 
depends strongly on it? 

The answer to this question is illuminated by considering 
the Usadel equations in the ferromagnetic region in limiting 
cases. From Eq. ([T5l) , we obtain several important properties: 

• In the presence of uniaxial spin-flip scattering (Sz = 
3, Sxy = 0), by taking the limit of r^f ^ we get 

— — ^- Then both singlet and triplet components 
are suppressed. This is also the case for isotropic spin- 
flip scattering (Sz = Sxy = 1). 

• With in-plane spin-flip scattering (Sz = 0), we get 0^ = 
Ol when taking the limit Tgf 0. Then, the singlet 
component should vanish as seen from definition in Eq. 
([271). 

• In the presence of pure spin-orbit scattering (Sz = 
Sxy = 0), by taking limit of Tso ^ we get 0^ = —O^. 
Then, the triplet component should vanish as seen from 
the definition in Eq. dZTl) . 

Therefore, we find that uniaxial and isotropic spin-flip scat- 
tering is harmful to both singlet and triplet componets while 
in-plane spin-flip scattering and spin-orbit scattering are detri- 
mental to the singlet and triplet components, respectively. 

These properties are also obtained by using linearized Us- 
adel equations in the ferromagnetic region, which may be for- 
mally obtained from Eq. ilSi by letting 6a- fa with a = ±, 
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FIG. 6: (Color online). Plot of the DOS of a ferromagnet/superconductor junction for a fixed exchange field /i/A = 0.3 which gives a ZEP in 
the absence of spin-dependent scattering. We fix d/^ = 1.0 and r = 0.5. 
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and assuming that |/± | <C 1: 



5± = 



1 r 
1 r 



2\{e ±h) - 



g^fSz 



2i{e ± h) 



g^f(2Sxy + Sz) 



(29) 



(28) 



As seen, the spin-orbit scattering rate only enters in the coef- 
ficient associated with the triplet component. Notice that in 
absence of an exchange field which renders ft = 0, the 
linearized Usadel equation is independent of the spin-orbit 
scattering rate. On a microscopical level, it is clear that the 
different dependence on spin-orbit scattering for the singlet 
and triplet anomalous Green's functions originates from the 
fundamental symmetries of these wavef unctions. It should be 
emphasized that both the singlet and triplet components are 
strongly affected by magnetic impurities, i.e. spin-flip scatter- 
ing, as seen from Eq. (|29l ). 

The above analysis emphasizes the importance of distin- 
guishing between different types of spin-dependent scattering 
in terms of understanding the behaviour of the DOS in a ferro- 
magnet/superconductor bilayer. In particular, we have shown 
that the effect of spin-orbit scattering may differ fundamen- 
tally from spin-flip scattering (originating e.g. from magnetic 
impurities), and that this is manifested in the interplay be- 
tween the singlet and triplet anomalous Green's function in 
the ferromagnetic layer. 



6. Minigap and spin-orbit scattering 

Finally, we consider the effect of spin-dependent scattering 
on the interesting manifestation of the proximity-induced su- 
perconducting correlations in a normal/superconductor junc- 
tion: the minigap^^. We first briefly recapitulate the results 
in the paramagnetic case. When a normal metal is placed in 
close proximity to a superconductor, a minigap opens up 
in the DOS of the normal metal. This minigap roughly scales 
like ~ Stt/'^r, and was originally studied by McMillan 
in a tunneling model of a normal/superconductor junction^. 



12 



Spin-orbit scattering 




0.4 0.6 0.8 



FIG. 8: (Color online). Plot of the singlet (lm{fs}) and triplet 
(Re{/t}) anomalous Green's functions induced in the ferromagnet 
for £ = right at the N/F interface (x = 0). We have used r = 0.5 
and d/^ = 1.0. 



The minigap is defined as an almost complete suppression of 
the quasiparticle DOS in a given energy interval. The minigap 
is a consequence of the effective backscattering that quasipar- 
ticles incident on the superconducting interface from the nor- 
mal side experience due to the presence of impurities. Con- 



sequently, the probability for transmission increases such that 
the DOS is nonzero for < e < A. For nearly perfect 
transparency of the interface, A^ can become close to A in 
magnitude. For a tunneling barrier r <C 1, the minigap be- 
comes very small. 

An interesting issue is how the minigap is affected by spin- 
dependent scattering. We provide numerical results to elu- 
cidate this question in Fig. [9l As seen from the plots, the 
influence of uniaxial spin-flip scattering is a gradual suppres- 
sion of the proximity-induced features in the DOS, consistent 
with previous results^^'^^'^^'^^'^^. The effect of isotropic spin- 
flip scattering is virtually identical to the uniaxial case. How- 
ever, quite surprisingly, the minigap shows strong resiliance 
towards increasing spin-orbit scattering as seen in Fig. [9j In 
fact, we find that even upon increasing g to values A, the 
minigap persists in the DOS as long as the exchange field is 
absent. Upon increasing the exchange field to values h A, 
the minigap slowly begins to close when the spin-orbit scat- 
tering rate becomes large. This feature may be understood 
by again resorting to Fig. [7] and Fig. [8l where a detailed 
study of the singlet and triplet components of the anomalous 
Green's function was conducted. In general, the spin-orbit 
scattering was shown to be detrimental for the triplet compo- 
nent while the singlet component still remained at a consider- 
able magnitude even for ^ 1. However, the triplet compo- 
nent was strengthened upon increasing the exchange field h, 
which explains why the minigap would close for larger val- 
ues of h upon increasing the scattering rate. This behaviour 
should be contrasted with spin-flip scattering, whether it is 
isotropic or uniaxial, which invariably suppresses both the sin- 
glet and triplet components of the proximity-induced anoma- 
lous Green's function. 
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FIG. 9: (Color online). Plot of the DOS at x = for a normal metal/superconductor junction using d/^ = 1.0 and r = 0.5. As seen, spin-flip 
scattering closes the minigap, whether it is uniaxial or isotropic. However, the minigap shows strong resiliance against increasing spin-orbit 
scattering. 



B. Josephson current lowing, fix Rb / Rf = 5 where Rb is the resistance of the 

Having concluded our study of the DOS in the F/S bilayer, 
we now turn to the supercurrent in a S/F/S junction. In the fol- 
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barrier at the interfaces and Rf is the resistance of the diffu- 
sive ferromagnetic layer. Also, we consider an exchange field 
h/A = 10, which should correspond to a typical weak ferro- 
magnetic alloy such as Cui-a^Nia^. In the following, we are 
particularly interested in examining the possibility of obtain- 
ing a strong deviation from the usual sinusoidal current-phase 
relationship. This may lead to the opportunity of creating a 
supercurrent- switching device, which we elaborate on below. 
We remind the reader that all the quantities used below were 
introduced and defined in Sec. HIl 

Consider first dependence of the /ci^ at -product on the nor- 
malized junction width d/^, where ^ = y^D/ 27rTc is the su- 
perconducting coherence length. We contrast the results ob- 
tained with a barrier of intermediate transparency {Z = 0.5) 
and low transparency {Z = 5.0) as the experimentally most 
relevant cases, shown in Fig. [TOl As seen, the intermediate 
transparency allows for a finite residual value of the supercur- 
rent at the O-TT transition point, where the first harmonic of 
the current-phase relationship vanishes. The effect of uniax- 
ial spin-flip scattering is seen to be a reduction of the residual 
value of the supercurrent. Even for intermediate transparency 
of the barrier, the residual current is reduced to immeasurable 
values for ^ ~ 5, where ^ = 7/A is a measure of the uniaxial 
spin-flip scattering. This issue has not been addressed previ- 
ously in the literature, and is useful to complement previous 
qualitative predictions with a more realistic quantitative anal- 
ysis. Also, it is seen from Fig. [TOlthat the transition points 
are translated towards higher junction widths upon increasing 
9- 

We next consider the temperature-dependence of the IqRn- 
product in Fig. [TTl In the wide-junction regime £ = St / A = 
0.05, the magnitude of IcRn is quite small and a O-tt tran- 
sition occurs for zero spin-flip scattering at T/Tc ~ 0.55. 
However, increasing the value of g is seen to vanish the O-tt 
transition completely. This behaviour is very distinct from the 
(i/^-dependence shown in Fig. [Tol where increasing spin-flip 
scattering reduces the residual value of the supercurrent at the 
transition point, but does not remove the transition completely. 
Consider now a shorter junction here modelled by f = 9.0, 
corresponding to (i/^ 0.6, for an intermediate barrier trans- 
parency Z = 0.5. Close examination of the transition point 
where the current changes sign reveals a small but finite resid- 
ual value of the supercurrent, in this case for T/Tc ~ 0.65. At 
first sight, the effect of increasing g then appears to amount to 
a complete removal of the O-tt transition point, rather than a 
suppression of the residual value as in Fig. [TOl However, 
by increasing g in smaller steps as shown in Fig. fT2h ) with 
S = 10.0, it is seen that the O-tt transition gradually van- 
ishes. This means that the transition point when considering 
the temperature-dependence is much more sensitive to spin- 
flip scattering than the width-dependence of Fig. \TU\ It should 
be noted that the O-tt oscillations vanish upon increasing g for 
the particular choice of the width d (corresponding to a cer- 
tain £) used here. For another choice of d, one might expect 
to have O-tt introduced upon increasing g. The main point 
is nevertheless that the temperature-dependence of the criti- 
cal current oscillations appear to be more sensitive to spin-flip 
scattering than the width-dependence of the same oscillations. 




FIG. 10: (Color online) Plot of junction- width dependence of the 
/ci? AT -product for an S/F/S Josephson junction in the intermediate 
{Z = 0.5) and low {Z = 5.0) barrier transparency regime. The 
temperature is fixed at T/Tc = 0.1. For each case, the effect of 
increasing uniaxial spin-flip scattering is seen to severly reduce the 
residual value of the critical current at the O-tt transition points. 



The origin of a residual value of the current at the sign- 
reversal point is a deviation from a purely sinusoidal current- 
phase relationship. To illustrate the strong deviation from a 
pure sinusoidal phase-dependence for the S = 9.0 case, con- 
sider Fig. [T2b). The O-tt transition can clearly be discerned 
from the plot. Right before the transition (T/Tc = 0.6), the 
maximum value of the current occurs at a negative value for 
IqRn- After the transition (T/Tc = 0.7), the maximum value 
occurs for a positive value of IqRn- For comparison, we have 
plotted a pure sinusoidal phase-dependence. Fig. fT2b ) shows 
that that higher harmonics [in this case sin(2(/))] dominate near 
the sign reversal point, giving rise to the residual value of the 
current. 
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FIG. 11: (Color online) Plot of the temperature-dependence of the 
/ci? AT -product for an S/F/S Josephson junction in the intermediate 
{Z = 0.5) and low {Z = 5.0) barrier transparency regime for dif- 
ferent values of the normalized Thouless energy 8. 



IV. DISCUSSION 

In our calculations, we have neglected the depletion of the 
superconducting and magnetic order parameter in the vicin- 
ity of the F/S interface. In the low transparency case, this 
is permissable^^. In the high transparency case, the depletion 
may be substantial. However, the qualitative features obtained 
in the DOS are known to be the same even if this depletion is 
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FIG. 12: (Color online) a) Plot of the temperature-dependence of 
the /ci? AT -product for weakly increasing spin-flip scattering, b) Plot 
of the current-phase relationship for the /ci? at -product for an S/F/S 
Josephson junction, illustrating the strong deviation from the sinu- 
soidal dependence and the O-tt transition. 



taken into account - the characteristic features are often sim- 
ply shifted in energy from the bulk value e = A to a reduced 
value A < A. Our approximation consisting of using the 
bulk Green's function in the superconducting regions is well 
justified also if we assume that the superconducting region is 
much less disordered than the ferromagnetic layer. Also, we 
have employed an effectively one-dimensional model to ac- 
count for the proximity-effect. In the present system, we are 
concerned with isotropic order parameters (s-wave supercon- 
ductivity), such that performing the same calculations in 2 or 3 
dimensional will yield similar results since the averaging over 
angles will have virtually no effect. 

Let us also comment on our choices of Thouless energy in 
the ferromagnetic region. The actual thicknesses employed in 
experiments on both F/S bilayers and S/F/S Josephson junc- 
tions vary greatly, but may in some cases be as small as < 20 
j^j^i9,20 Nq^^ typical superconducting coherence length 
also varies a lot, ranging from 38 nm in Nb to 1600 nm in Al. 
Since the coherence length is defined as ^ = D / /S. with the 
diffusion constant D in the ferromagnetic region, we find that 
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if £t/ A = X, then the thickness d of the ferromagnetic region 
is given as d = ^/ y/x. Assuming a superconducting mate- 
rial with ^ = 100 nm, choosing e^/A = 0.1 corresponds to 
(i :^ 315 nm. In order to make contact with the experimen- 
tal situations that employ a ferromagnetic layer thickness of 
d 20 nm, we would then need a value of x = 25. This 
is the motivation for our choice of the large Thouless energy 
St/^ = 100, corresponding to d/^ = 0.1. 

Finally, we discuss the possible realization of a 
supercurrent- switch. The prospect of obtaining a dissi- 
pationless current- switching device relies on the opportunity 
to achieve sufficiently large residual values of the supercur- 
rent at the O-TT transition points. As seen from our results, 
this may be obtained for the intermediate barrier transparency 
regime, which agrees with the qualitative conclusions of 
Ref. [Isl. However, the effect of uniaxial spin-flip scattering 
has a detrimental effect on this residual value. While the 
width-dependence (Fig. (TO]) shows some resilience towards 
an increased concentration of magnetic impurities, the 
O-TT transition is highly sensitive to spin-flip scattering in 
temperature-dependence (Fig. [TT]). To switch the sign of the 
current, one has to increase either the width or temperature of 
the junction by an infinitesimal amount right at the transition 
point. Since only temperature can be manipulated in this way 
in a realisitic experiment, we arrive at the conclusion that 
the realization of a current- switch device in diffusive S/F/S 
junctions relies on samples with high quality interface and 
very low amounts of magnetic impurities. 

V. SUMMARY 

In conclusion, we have numerically studied the local den- 
sity of states in a proximity ferromagnet/superconductor 
structure in the dirty limit. Our results take into account an ar- 
bitrary rate of spin-flip scattering and arbitrary interface trans- 
parency, focusing on a moderately transparent interface. This 
regime cannot be reached in the standard limiting cases of an 
ideal or tunneling interface. We have studied three cases for 
the size of the exchange field h compared to the supercon- 
ducting gap A: i) A, ii) h > A, in) h:^A.ln each case, 
we considered several values of the thickness of the ferromag- 
netic layer, barrier transparencies, and also different types of 
spin-dependent scattering to obtain the energy-resolved DOS. 
In doing so, we have clarified characteristic features that may 
be expected in the various parameter regimes accessible for 
the ferromagnetic film. Since our results take into account ar- 
bitrary proximity effect and magnetic impurity concentration, 
they should serve as a useful tool for a quantitative analysis 
of experimental data. In particular, we investigated the effect 
of spin-dependent scattering on the zero-energy behaviour ob- 
served in the DOS, which displayed the full range from a fully 
developed minigap to a peak- structure. By analyzing in detail 
the singlet and triplet part of the anomalous Green's function, 
we come to the important conclusion that it is necessary to 
distinguish between different types of spin-dependent scatter- 
ing in order to correctly interpret DOS-measurements in fer- 



romagnet/superconductor bilayers. Specifically, we find that 
the effect of spin-flip scattering (both uniaxial and isotropic) 
may differ fundamentally from the effect of spin-orbit scatter- 
ing. The reason for this is that for weak exchange fields, the 
singlet component of the anomalous Green's function remains 
essentially unaltered by spin-orbit scattering, while the triplet 
component is strongly suppressed. 

We have also investigated the supercurrent in diffusive su- 
perconductor/ferromagnet/superconductor junctions, allow- 
ing for arbitrary concentration of magnetic impurities and ar- 
bitrary interface transparency. We have investigated the effect 
of spin-flip scattering on the residual value of the supercurrent 
at the O-TT transition points, and find a much weaker sensitivity 
to magnetic impurities in the width-dependence compared to 
the temperature-dependence of the at -product. We have 
proposed that a finite and measurable residual value of the su- 
percurrent may be obtained in the intermediate barrier trans- 
parency regime, although spin-flip scattering has a detrimen- 
tal effect on this residual value. For samples with high qual- 
ity interface and very low concentration of magnetic impu- 
rities, the residual value may be exploited to obtain efficient 
supercurrent- switching simply by altering the temperature of 
the system. 
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APPENDIX 

We here define the star-product which enters the Eilenberger 
equation Eq. For any two functions A and B, we have 

where the differentiation operators denote derivation with re- 
spect to the variables T and e in the mixed representation. 
Note that if there is no explicit time-dependence in the prob- 
lem, the star-product reduces to regular multiplication. 
The Pauli-matrices used in this paper are defined as^ 

a= (! J)' Z2 = (j o')' ^= (o -i)' 

?)'^=(Qi)'^^ = (t i)' 
^1 = (a t) ' = (4 ~Q~) ' = (q "i) ■ 

(31) 
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